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Two fundamental facts of the modern wave turbulence theory are 1) existence of power energy 
spectra in fc-space, and 2) existence of "gaps" in this spectra corresponding to the resonance cluster- 
ing. Accordingly, three wave turbulent regimes are singled out: kinetic, described by wave kinetic 
equations and power energy spectra; discrete, characterized by resonance clustering; and mesoscopic, 
where both types of wave field time evolution coexist. In this paper we study integrable dynamics 
of resonance clusters appearing in discrete and mesoscopic wave turbulent regimes. Using a novel 
method based on the notion of dynamical invariant we establish that some of the frequently met 
clusters are integrable in quadratures for arbitrary initial conditions and some others - only for par- 
ticular initial conditions. We also identify chaotic behaviour in some cases. Physical implications of 
the results obtained are discussed. 
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I. INTRODUCTION 

The broad structure of modern nonlinear science born 
at the edge of physics and mathematics includes an enor- 
mous number of applications in cosmology, biochemistry, 
electronics, optics, hydrodynamics, economics, neuro- 
science, etc. The emergence of nonlinear science itself as 
a collective interdisciplinary activity is due to the aware- 
ness that its dynamic concepts first observed and un- 
derstood in one field (for example, population biology, 
flame-front propagation, non-linear optics or planetary 
motion) could be useful in others (such as in chemical 
dynamics, ncuroscience, plasma confinement or weather 
prediction). The theory of integrable Hamiltonian sys- 
tems, a generalization of the classical theory of differ- 
ential equations, is the nucleus of the whole nonlinear 
science. Various classifications of integrable systems are 
presently known which turned out to be quite useful for 
physical applications. Classifications are known based on 
the various intrinsic properties of integrable systems [53] : 
symmetries, conservation laws, Lax-pairs, etc. In [3] the 
general classification of integrable Hamiltonian systems 
is presented based on the form of their topological invari- 
ants. The usefulness of this classification is demonstrated 
in several problems on solid mechanics. In particular, it 
is proven that two famous problems - the Euler case in 
rigid body dynamics and the Jacobi problem of geodesies 
on the ellipsoid- are orbitally equivalent. In [12] the idea 
of classification is presented based on normal forms of 
a certain class of bi-hamiltonian PDEs. Miscellaneous 
hierarchies of integrable PDEs are presented in [48] . 

The list can be prolonged further but the main point 
for us presently is the following: the notion of integrabil- 
ity itself is ambitious! There are many quite different def- 
initions of integrability, for instance intcgrability in terms 
of elementary functions (equation y = —y has the explicit 
solution y = a sin (a; + b) ); integrability modulo class of 
functions (equation y = f(y) has general solutions in 
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FIG. 1: Color online. Schematic representation of wave turbulent regimes. 



terms of elliptic functions), etc. An example of less obvi- 
ous definition of integrability is C-integrability, first in- 
troduced in [6]: integrability modulo change of variables, 
meaning that a nonlinear equation is called C-intcgrable 
if it can be turned into a linear equation by an appropri- 
ate invertible change of variables. For instance, Thomas 
equation tp xy + aip x + /3ip y + ip x ip y = is C-integrable. 
Profound discussion on the subject can be found in [34]. 
In the present paper, integrability is interpreted in terms 
of the existence of a number of independent dynamical 
invariants of the system; for each in-this-sense-integrablc 
system, solutions are then written out in quadratures. 

The dynamical systems we are interested in, describe 
nonlinear resonance clusters appearing in evolutionary 
dispersive PDEs in two space variables. Nonlinear reso- 
nances are ubiquitous in physics. They appear in a great 
amount of typical mechanical systems [13, 41], in engi- 
neering [8, 17, 42, 57], astronomy [40], biology [15], etc. 
etc. Euler equations, regarded with various boundary 
conditions and specific values of some parameters, de- 
scribe an enormous number of nonlinear dispersive wave 
systems (capillary waves, surface water waves, atmo- 
spheric planetary waves, drift waves in plasma, etc.) all 
possessing nonlinear resonances. 

The classical approach of statistical wave turbulence 
theory in a nonlinear wave system assumes weak non- 
linearity, randomness of phases, infinite-box limit, and 
existence of an inertial interval in wavenumber space 
(fco,fci), where energy input and dissipation are bal- 
anced. Under these assumptions, the wave system is 
energy conserving, and wave kinetic equations describ- 
ing the wave spectrum have stationary solutions in the 
form of Kolmogorov-Zakharov (KZ) energy power spec- 
tra k~ a , a > 0, ([55, 63, 66], etc.). 

As it was first established in the frame of the model 
of laminated turbulence, [22], KZ-spectra have "gaps" 
formed by exact and quasi-resonances (that is, resonances 
with small enough resonance broadening). This yields 
two distinct layers of turbulence in an arbitrary nonlin- 
ear wave system - continuous and discrete - and their in- 



terplay generates three possible wave turbulent regimes: 
kinetic, discrete and mesoscopic as it is shown in Fig. 1 
(see [10] for more discussion). 

The existence of mesoscopic regime has been first con- 
firmed in numerical simulations with dynamical equa- 
tions for surface gravity waves in [67], while discrete 
regime has been first described in [25]. These theoret- 
ical findings are confirmed by numerous laboratory ex- 
periments. For instance, in the experiments with gravity 
surface wave turbulence in a laboratory flume, [11], only 
a discrete regime has been identified while in [62] coex- 
istence of both types of time evolution has been estab- 
lished. Taking into account additional physical param- 
eters in a wave system transition from kinetic to meso- 
scopic regime can be observed as it was demonstrated in 
[9] for capillary water waves, with and without rotation. 

From a mathematical point of view, the very special 
role of resonant solutions has been first demonstrated by 
Poincare who proved, using Calogero's terminology, that 
a nonlinear ODE is C-integrablc if it has no resonance so- 
lutions (see [2] and refs. therein). This statement allows 
the following Hamiltonian formulation [63]: 

ictk = dH/da^, (1) 

where Ok is the amplitude of the Fourier mode corre- 
sponding to the wavevector k and the Hamiltonian T-L is 
represented as an expansion in powers 7-Lj which are pro- 
portional to the product of j amplitudes a^. Then the 
cubic Hamiltonian has the form 

H3 = V 1 3 2 aia2a3^23 + complex conj., 

ki,k 2 ,k 3 

where for brevity we introduced the notation aj = 
and 623 = 5(ks — ki — k2) is the Kronecker symbol. 
If H3 7^ 0, three-wave resonant processes are dominant. 
These satisfy the resonance conditions: 

w(k x ) + cj(k 2 ) - w(k 3 ) = 
ki + k 2 - k 3 = 0, U 
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where w(k) is a dispersion relation for the linear wave 
frequency. Further on, the notation u)k is used for w(k). 
The corresponding dynamical system has a general form 

iBk = (V^BiB^ + 2V^*B 1 B;S 1 k3 ) (3) 

ki,k 2 

(notations Bj are used further on for the resonant 
modes). If = 0, four- wave resonances have to be 
studied, and so on. To confirm that =^ and three- 
wave resonances are dominant, one has to find solutions 
of (2) and check that V^ 2 ^ at least at some resonant 
triads. Afterwards the corresponding dynamical system 
has to be studied. 

It has been first proven in [21] that for a big class 
of physically relevant dispersion functions uj, the set of 
all wavevectors satisfying (2) can be divided into non- 
intersecting classes and solutions of (2) can be looked 
for in each class separately. The method of g-class de- 
composition first introduced in [23] has been developed 
specially for solving systems of the form (2) in integers; 
details of its implementation for various rational and irra- 
tional dispersion functions are given in [27] -[29]. General 
description of the q-class method and corresponding pro- 
gramming codes are given in [26], in Ch.3 and Appendix 
correspondingly. 

An immediate consequence of the g-class method is 
that dynamical system (3) can be reduced to a few 
dynamical systems of smaller order, and each of these 
smaller dynamical systems can be investigated indepen- 
dently from all others. In [32], construction of a set of re- 
duced dynamical systems corresponding to the solutions 
of (2) and the systems themselves are given explicitly (as 
an example, resonances of oceanic planetary waves were 
considered in the spectral domain < m,n < 50). The 
integrability of some resonance clusters has been studied 
in [4]. 

The main goal of the present paper is to study system- 
atically the integrable dynamics of the most frequently 
met resonance clusters. We begin with a brief introduc- 
tion of NR-diagrams (NR for nonlinear resonance) which 
give a handy graphical representation of a generic reso- 
nance cluster and allow us to recover uniquely the dy- 
namical system corresponding to each cluster [26]. 



II. NR-DIAGRAMS 

In systems with cubic Hamiltonian, a resonant triad is 
called primary cluster [29] (a resonant quartet is a pri- 
mary cluster in a system with quadric Hamiltonian, and 
so on) . All other clusters (formed by a few primary clus- 
ters connected via one of a few joint modes) are called 
generic clusters or simply clusters. The dynamical sys- 
tem for a complex triad in the standard Manley-Rowc 
form (that is, with one interaction coefficient Z) reads 

B\ = ZB 2 B 3 , B 2 = ZB* x B 3l £> 3 = — ZB\B 2 , (4) 



as) j m /yWyp" 




FIG. 2: Topological structure of the cluster set for the oceanic 
planetary waves, lo ~ l/\/m 2 + n 2 , in the domain m, n < 50. 
7 types of resonance clusters have been found, the number of 
the clusters of each type is shown in parenthesis. 

and is known to be integrable (e.g. [61]), with two con- 
servation laws in the Manley-Rowe form being 

ha = \B 2 \ 2 + \B 3 \ 2 , J 13 = l^xl 2 + |B 3 | 2 . 

Due to the criterion of nonlinear instability for a triad 
[16], the mode with maximal frequency, U3, is unstable 
while the modes u\ and u 2 are neutral. Originally, this 
fact has been deduced directly from the equations of mo- 
tion, [16], but it can easily be seen from the the form of 
Manlcy-Rowe constants, [31]. 

This means that the form of dynamical systems and 
accordingly time evolution of the modes belonging to a 
generic cluster depends crucially on the fact whether joint 
modes within a cluster are stable or unstable. With the 
purpose to distinguish between these cases, the notations 
A-mode (active) and P-mode (passive) are introduced for 
w 3 -mode and w\- and W2-modes respectively, [31]. This 
allows to describe all possible connection types within a 
generic cluster. For instance, 1-mode connection of two 
triads can be of AA-, AP- and PP-type; 1-mode connec- 
tion of three triads can be of AAA-, AAP-, APP-type and 
PPP-type, 2-mode connection between two triads can be 
of AA-PP-, AP-AP-, AP-PP- and PP-PP-types, and so 
on. 

In the topological representation [32] of the solution 
set of (2) this dynamical information has been kept un- 
cxplicit, as part of a programming code used to construct 
dynamical system, while each triad within a cluster was 
shown as an unmarked triangle (see Fig. 2). In this rep- 
resentation each vertex shown as a circle denotes one 
resonant mode. This representation has been slightly im- 
proved in [31] where each ^4-mode has been marked by 
two arrow-edges coming from ^4-mode to both P-modcs 
in each triad. However, for a larger clusters this repre- 
sentation becomes too nebulous. 

More compact graphical representation of a resonance 
cluster is given by its NR-diagram first introduced in [25] , 
both for three- and four- wave resonance systems. In a 
NR-diagram each vertex represents not a resonant mode 
but a primary cluster, that is, a triad and a quartet in a 
three- and four-wave system correspondingly 
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FIG. 3: NR-diagrams for some resonance clusters shown in 
Fig.2. 



A NR-diagram in systems with cubic Hamiltonian con- 
sists of following building elements - a triangle and two 
types of half-edges, bold for A-modc and dotted for P- 
mode. It can be proven ([26], Ch.3) that in this case 
the form of NR-diagram defines uniquely corresponding 
dynamical system. Examples of NR-diagrams for some 
resonance clusters shown in Fig.2 are displayed in Fig. 3. 

Below, dynamical systems are given for three generic 
clusters shown in Fig.2, on the left (dynamical systems 
for the other clusters are omitted for sake of space) : 

1. Cluster consisting of two triads a and b, whose con- 
necting mode is active in one triad and passive in the 
other triad, say B 3a = Bib- In other words, a cluster 
with one AP-connection. It is called AP-butterfly [31] 
and its dynamical system is 

B\ a = Z a B 2a B 3a , B 3 b = —ZtB 3a B 2 b , 
< B 2a = Z a B* a B 3a , B 2 b = Z b B^ a B 3b , (5) 
B 3a = — Z a B\ a B 2a + ZbB^Bst, 

2. Cluster consisting of three triads a, b and c, with 
one AA- and one PP-connections, say, B 3a = B 3 b and 
Bib = Bic- The dynamical system reads 





= Z a B2 a B 3a , 


B 2a = Z a B* a B 3a 


B 3a 


= —ZuB\aB<ia 


— ZbB\bB 2 b 


Bib 


= ZbB^B^a + 


ZB 2c B 3c , 


B 2b 


= ZbB\ h B 3ai 


B 2c = Z c Bl b B 3c , 


,B 3c 


= — Z c B\bB 2c . 





3. Cluster consisting of four triads a, b, c and d, with 
two AA- and one PP-connections, say, B 3a = B 3 b = B 3c 
and B\ c = Bid- The corresponding dynamical system is 



of the form 

&\a — 6 a D 2a 3a , D 2a — /j a & la U 3a} 

B 3a = —Z a B la B 2a — Z c B lc B 2c — ZdBidB 2c i, 
^ B\b = ZbB2 b B 3a , B 2 b = ZbB\ b B 3a , 
B 3 b = —ZbB\bB 2 b 
Bic = Z c B 2c B 3a + ZdB2 d B 3 d, 
K B 2c = Z c B{ c B 3a , B 2 d = ZdB\ c B 3 d- 

The Manley-Rowe constants can be written out imme- 
diately for each of these systems, being combinations of 
corresponding constants for each triad. For instance, for 
(5) they have the form 

{-^126 = |-Bih| 2 — |i?2fc| 2 , h3b = \B 2 b\ 2 + \B 3b \ 2 , , . 
Iab = \B lb \ 2 + \B 3a \ 2 + \B 3b \ 2 . U 

NR-diagrams also give us immediate qualitative infor- 
mation about the energy percolation within bigger res- 
onance clusters, for the case when two P-modes, form- 
ing a PP-connection, have small initial amplitudes com- 
pared to the amplitudes of A-modes in the connected 
triads. It follows then from the Hasselmann's criterion 
of instability [16] that a PP-connection can be regarded 
as an obstacle for the energy percolation in both direc- 
tions and this cluster can be then regarded practically as 
two independent triads. Analogous considerations show 
that AA-connection allows energy percolation in two di- 
rections while AP-connection - in one direction. In this 
sense one can define, for certain initial conditions, PP- 
rcductions of resonance clusters, whereby a cluster is ap- 
proximated by smaller clusters obtained by cutting off 
the PP-connections from the original cluster. 

While planning laboratory experiments, one has to be 
very careful with these theoretical findings. For instance, 
in the experiments reported in [7] , a chain- like cluster of 
three connected triads has been identified with one PP- 
connection. However, this connection could not be disre- 
garded while a very small "parasite" frequency generated 
by electronic equipment was enough for initiating the en- 
ergy exchange among the modes of all three triads (sec 
[26], Ch.4, for detailed explanations). 

The main difference between NR-diagram and statis- 
tical diagrams used in wave turbulence theory and origi- 
nated from Feynman diagrams can be formulated as fol- 
lows. Each statistical diagram corresponds to one term 
in the asymptotic expansion and does not allow to com- 
pute the amplitudes of the scattering process. On the 
other hand, a NR-diagram describes completely a reso- 
nance cluster and allows to write out explicit form of the 
dynamical system on the modes' amplitudes. 

As it will be shown below, connection types within a 
cluster define indeed the integrability of the correspond- 
ing dynamical systems. In order to demonstrate it we will 
use the notion of dynamical invariant first introduced in 
[4] which is given in the next section and illustrated by 
the example of harmonic oscillator. 
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III. DYNAMICAL INVARIANTS 
A. Definition 

From here on, general notations and terminology will 
follow Olver's book [53] and Einstein convention on re- 
peated indices and = df jdx % . Consider a general N- 
dimensional system of autonomous evolution equations 
of the form: 

^(t) = AV(t)), i = l,...,N. (9) 

Any scalar function f(x % ,t) that satisfies 

|(/(^(*),*)) = |/ + AV, = 

is called a conservation law in [53]. It is easy to see that 
this definition gives us two types of conservation laws: 
(i) those of the form f{x l ) (no explicit time-dependence), 
and (ii) those of the form f(x l ,t), where the time depen- 
dence is explicit. The first type determines an invariant 
manifold for the dynamical system (9) (time-independent 
conservation law) and the second type constrains the 
time evolution of the system within the invariant man- 
ifold^) (time- dependent conservation law). To keep in 
mind the difference between these two types of conserva- 
tion laws, we call the first type just a conservation law 
(CL), and the second type - a dynamical invariant. 

We are interested in determining the solution 
x l (t), i = 1, . . . , N, of a given dynamical system of 
the form (9). One possible way to do that is by find- 
ing N functionally independent dynamical invariants for 
the system (9). This is equivalent to finding (N— 1) func- 
tionally independent conservation laws and one dynami- 
cal invariant (the equivalence can be proven, for example, 
using the implicit function theorem). 

As it was shown in [5] , in some cases the knowledge of 
only (N — 2) functionally independent CLs is enough for 
constructing explicitly: (i) a new CL functionally inde- 
pendent of the others, and (ii) a corresponding dynamical 
invariant, determining the solution £*(£), i = 1, . . . , N. 
This follows from the Theorem on (N — 2)-integrability 
[5], whose formulation is given below for the readers' con- 
venience. 

Theorem on (N — 2)— integrability. Let us assume 
that the system (9) possesses a standard Liouville volume 
density 

ptf): ( P A% = 0, 

and (N — 2) functionally independent CLs, 
H 1 ,...,^- 2 . Then a new CL in quadratures can 
be constructed, which is functionally independent of the 
original ones, and therefore the system is integrable. 



B. Example: damped harmonic oscillator 

1. Dynamical invariants, CLs and solutions 

To illustrate the complementarity of conserved laws 
and dynamical invariants, we present an illustrative ex- 
ample from mechanics for the case N = 2. Consider the 
damped harmonic oscillator. The equations of motion in 
non-dimensional form can be written as: 

q=p, p = -q-ap, (10) 

where a > is the damping coefficient. This is a dynam- 
ical system of the form (9) with N = 2. 

Now we want to fully determine the solution of the 
dynamical system (10). For this we need to know both 
a CL and a dynamical invariant. Indeed, let us consider 
separately the cases a = (harmonic oscillator) and < 
a < 2 (sub-critically damped harmonic oscillator). 

1. Case a — 0. We have the CL 

E(q,p) = 1/2 (p 2 +q 2 ) (11) 
(energy) and the dynamical invariant 

T(q,p,t) = t - arctan(g/p). (12) 

Since 

j t (E(q(t),p(t))) = , l(T(q(t),p(t),t)) = 0, 

then we have 

E(q(t),p(t)) = E Q , T(q(t),p(t), t) = T , 

constants depending on the initial conditions q(0),p(0). 
This information is enough to find the solution q{t),p(t) 
of the system: 

q(t) ^Vmism(t-T ) 
p{t) = \/2£b cos(t - T ) 

which can be checked by direct substitution in (10). 

In Fig. 4, upper panel, we show level surface of conser- 
vation law E = 30 (on the left), level surface of dynamical 
invariant T = (middle) and solution trajectory q(t),p(t) 
(on the right). This solution is actually the intersection 
of the level surfaces E = 30 and T = 0; for completeness 
of presentation, we show together the level surfaces and 
the solution in Fig. 5 (left panel). 

Notice that coordinates q, p are not suitable for a global 
paramctrization of dynamical invariant T, (12), because 
arctan(x) is multi-valued. This problem can easily be 
overcome by using new coordinates (R,9): 

q = Rsin9, p = Rcos9 (14) 

which allow us to rewrite T as 

T(q,p,t)=f(R,9,t)=t-e, (15) 
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FIG. 4: Color online. Upper panel: harmonic oscillator (case a = 0): left upper panel, a level surface of conservation law 
E(q,p) — 30, Eq.(ll); middle upper panel, a level surface of dynamical invariant T(q,p,t) — 0, Eq.(15); right upper panel, 
solution trajectory (q(t),p(t)), Eq.(13), corresponding to E = 30, T = 0. Lower panel: sub-critically damped harmonic oscillator 
(case < a < 2): left lower panel, a level surface of conservation law C(q,p) — 30, Eq.(19); middle upper panel, a level surface 
of dynamical invariant D(q,p,t) = 20, Eq.(16); right upper panel, solution trajectory (q(t),p(t)), Eq.(20), corresponding to 
C = 30, D = 20. 




FIG. 5: Color online. Left panel: harmonic oscillator (case a = 0): combined plot of level surface of conservation law 
E(q,p), level surface of dynamical invariant T(q,p,t) and solution trajectory (q(t),p(t)). Notice the general property that the 
intersection of the level surfaces of E and T is the solution trajectory. Right panel: sub-critically damped harmonic oscillator 
(case < q < 2): combined plot of level surface of conservation law C(q,p), level surface of dynamical invariant D(q,p,t) and 
solution trajectory (q(t),p(t)). Notice the general property that the intersection of the level surfaces of C and D is the solution 
trajectory. 
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thereby eliminating the ambiguity. The plot of level sur- 
face T = in the middle upper panel of Fig. 4 was done 
using the parameters (R, t) . 

2. Case < a < 2. Let a/2 = simp, then a dynamical 
invariant for the system is known: 



coefficient a. A simple pendulum oscillating at small am- 
plitudes is easy to construct, consisting of a massive bob 
attached to the ceiling by a light (i.e. massless) string 
or rod. Damping can be easily introduced by attaching 
to the pendulum rod a sheet made of a light material. 
For small amplitudes, the angular frequency of the os- 



D(q,p,t) = {p 2 + q 2 + 2pqsinip) exp(2isiny>). (16) dilations is given by u = yJg/L, where g w 9.8[m/s 2 } is 



We need to find a CL in order to determine the solution. 
Here we simply state the following CL: 



C{q,p) 



cos ip 
2 

x cxp 



(P 2 



q 



2pq sin tp) 
M 



(17) 



2 tan tp arctan(— sec ip + tan tp) 
P 



The method of construction of this CL is not important 
right now, it will be detailed in the next section. In Fig. 4, 
lower panel, we show level surface of conservation law 
C = 30 (on the left), level surface of dynamical invariant 
D = 20 (middle) and solution trajectory q(t),p(t) (on the 
right). This solution corresponds to the intersection of 
the level surfaces C = 30 and D = 20; for completeness, 
the level surfaces and the solution are shown together in 
Fig. 5 (right panel). 

Similar to the case of a = dynamical invariant, the 
CL C(q,p) is not globally defined. The following change 
of variables is needed: 

q = R sec tp s'm(8 — ip), p = R cos 9 sec (p, (18) 

in order to parameterize globally the CL. The result is 

C{q,p) = C{R,9) = ^|^i? 2 exp(26»tan (/ 3). (19) 

It is important to realize that the two changes of variables 
(14) and (18) are suggested by the form of the respective 
invariants. Moreover, in the new variables (R,9), both 
T and C take a simpler form (see (15), (19)). The reason 
for it is clear in the case a = because (R, 9) are the 
well-known action-angle variables. In the general case, 
the variables (R, 9) determine a covering of the original 
variables. 

The solution of the dynamical system is finally 



q(t) = 
Pit) 



\f~D sin(t cos tp + h \ag{%) cot ip — tp) 



cos tp exp(i sin ip) 
D cos(i cos ip + | log(^) cot ip) 
cos tp exp(i sintyj) 



(20) 
(21) 



2. Description of a laboratory experiment 

To check the constancy of the conservation law C(q,p) 
in a school laboratory experiment, we give here detailed 
description of the possible experiment in the case of sub- 
critical damping < a < 2. Moreover, we provide evi- 
dence of the practical and physical interest of the conser- 
vation law C(q,p): it can be used to measure the damping 



the acceleration of gravity and L is the pendulum length. 
We assume ui to be known, though of course it could be 
measured experimentally. From now on, we choose units 
where uj = 1, so the damped pendulum equations for 
small amplitudes reduces to the system (10), where q is 
the position of the oscillating bob and p is its velocity. 

With the aid of cheap electronic detection equipment 
it is possible to measure accurately the time t, position q 
and velocity p of the oscillating bob at certain instances, 
namely when the bob passes near a detector. This data is 
sent to a computer database and analyzed with a software 
that comes along with the detection equipment. As a 
result one obtains a set of data points of the form 



(q n ,Pn,tn), n = 0,...,N 



(22) 



With the aid of only one detector, we can in principle 
measure data at instances when q n = 0. To wit, we first 
calibrate the detector by setting the q = coordinate as 
the equilibrium position of the mass. An experimental 
realization consists in producing a small amplitude os- 
cillation and acquiring data points of the form (22). By 
looking at (18), these instances correspond to 

9 n = tp + rm, p n = (-l) n R n , 

where tp is the unknown parameter related to the damp- 
ing coefficient by a = 2simp. From (19) we obtain data 
values 

cos^ o 

G„ = — — p n exp{2{tp + nir) tan tp). 

If C is a conservation law, then C n must be constant, 
independent of n. Therefore we should have 

p 2 n exp(2n7r tan tp) = const = p\. 

In practice one can measure tan tp by finding a linear fit 
of log \p n \ versus n: 

log|p„| = log |p 1 -rmt&nip, n=l,...,N. (23) 

Once tan tp has been obtained, the damping coefficient a 
is readily computed. 

Notice that one could use alternatively the dynamical 
invariant D(q,p,t) to compute a, from the data (22). 
Setting to = 0, from (16) we obtain 

D„ = p 2 n exp(2i„ sin<p) = const = p\. 

We can measure sin tp by finding a linear fit of log \p n \ 
versus t r . : 



log |Pn | = log|p | - t n smtp, n=l, 



,N. 



(24) 
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Finally, notice that these two ways to measure the damp- 
ing coefficient imply a third way: by comparing (23) and 
(24) we solve for the data times 

t n = nir I cos tp , n = 1, . . . , N. 

In this case, a simple linear fit of t n as a function of n will 
allow us to obtain tp and thereby the damping coefficient 
a. 

3. Construction of conservation law 

To illustrate the procedure of construction of a con- 
servation law, we take as an example the sub-critically 
damped harmonic oscillator, i.e. (10) with < a < 2. 
Here N = 2 and A 1 = p, A 2 = — q — ap. The dynami- 
cal system is just 2-dimensional and we will write it as a 
vector (A 1 , A 2 ) T . The Theorem requires the existence of 
a standard Liouville volume density p(x 1 ,x 2 ) satisfying 

(pA 1 ),! + ( P A 2 ), 2 = 0, (25) 

and does not require the knowledge of conservation laws. 
In general, a Liouville density, solution of (25), is inter- 
preted as follows. A small region lZ(t) with a volume 
V(t) in phase space (2:1,0:2), will evolve in time due to 
the dynamical system (10). Then, p is defined in such a 
way that the product pV(t) is conserved in time as lZ(t) 
evolves. For the harmonic oscillator, it is well known that 
the volume of 7Z(t) is preserved, i.e., a constant function 
is a Liouville density. For the damped harmonic oscil- 
lator, a direct check shows that a Liouville density is 

p(q,p) = (q 2 + p 2 + 2aqp) . With this information we 
just need to solve (26) for H: 

U ;)(£)-(#)■ (26 » 

The answer can be obtained by direct integration (see [4] 
for more details): 

fqa arctanf , 2p+qa )\ 

[q ' P> [ ^-(q 2 M + a 2 )) ) 

| \og(p 2 + q 2 +pqa) ^ 

The conservation law C(q,p) given by Eq.(17) is a func- 
tion of H , chosen for its nice form: 

C=^ ex p ( 2H). 

IV. TRIAD 

As it was shown above, the notion of dynamical invari- 
ant is an important tool for constructing new physically 



relevant conservation laws that can afterwards be stud- 
ied in a simple laboratory experiment. In this section we 
would like to use this approach to prove integrability of 
a complex triad with dynamical system (4). Of course, 
integrability of (4) is a well-known fact (e.g. [61]). How- 
ever, the explicit solution of (4) is usually written for 
a particular case, namely, when the dynamical phase - 
a phase combination corresponding to the chosen reso- 
nance conditions - is either zero or constant ([44], p. 132, 
Eq.(6.7); [54], p.156, Eq.(3. 26.19), etc.). 

On the other hand, it is well known that dynamical 
phases play a substantial role in the dynamics of reso- 
nant clusters, e.g. [58], and their effect can easily be 
observed in numerical simulations, [5]. This was our mo- 
tivation for constructing first an explicit solution in the 
amplitude-phase presentation, for (28). Thus, (4) is used 
for a preliminary check of our method. 

Another important point is the following. As it was 
shown in the papers [36, 37, 39], an elastic pendulum 
with suitably chosen parameters can be used as a me- 
chanical model of a resonant triad, and the results can 
be applied for the description of large-scale motions in 
the Earth's atmosphere. In fact, this simple mechanical 
model can be used for a laboratory study of dynamical 
characteristics of primary clusters in an arbitrary system 
with cubic Hamiltonian. This is why the explicit ana- 
lytical formulas for all dynamically relevant parameters, 
given below, are important. 



A. Integrability 

In this case the system can be reduced to N = 4 (see [4] 
for more details), the Theorem on (N — 2) -integrability 
can be applied and we obtain the following CL: 

H T = IMB1B2BI) , 

which is the canonical Hamiltonian for this case and can, 
of course, be written out directly. A dynamical invariant 
for this system was originally presented in [4], in terms 
of the three real roots R\ < R2 < R3 of the cubic poly- 
nomial 

x 3 +x 2 = 2/27 - (27Ht— 

(/l3 + /23)(/l3-2/2 3 )(/23-2/l3))/27(/ 2 3-/ 13 / 2 3 + / 2 2 3) 3/2 , 

but these roots' dependence on the coordinates or the 
CLs was not made explicit. Moreover, the explicit 
solution for the amplitudes Cj and phases 9j in the 
amplitude- phase representation Bj = Cj exp (iOj) was 
not provided. Here we improve the form of dynamical in- 
variant and also produce explicit and useful expressions 
for the full solution, based on the trigonometric repre- 
sentation of the three real roots in the so-called Casus 
Irreducibilis. 
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B. Amplitude-phase representation 



C. Solutions for amplitudes 



Sys.(4) in the standard amplitude-phase representation 
Bj = Cj exp(i6j) reads: 



C\ = ZC2C3 cos<£, 

C-l = ZC\C 3 COS<y9, 
C3 = —ZClC2 C0S(y9, 



(28) 



ip = -ZH T (Ci 



Co 



where ip — 9\ + 62 — 3 is the dynamical phase. The 
conservation laws (5) do not change their form in the 
new variables: 

I23 = Cf + Cf , I13 = C 1 + C 3 , 



but the Hamiltonian Ht reads now 

Ht = C1C2C3 siny;. 
Let us introduce new variables: 

P = hz/h-3 
and a € [0,7r] defined by 

(-2 + 3 p + 3p 2 -2p 3 ) I 13 3 - 27 H T 2 



cos a 



2(1 -p + f ?)*l 13 z 



(29) 



(30) 



(31) 



Notice that |cosa| < 1 for dynamically accessible sys- 
tem's configurations. Indeed, the use of intermediate 

3 

variables p = 2 (l — p + p 2 ) 2 and q = —2 + 3p + 3p 2 — 
2 p 3 , allows one to conclude immediately that p > and 
V > M V p . Both inequalities become equalities if p = 
or p = 1. This yields 



and 



(cosa) r 



(p + > , for cosa = — 1, 



q/p<l, for £T T = 0, 



where (H T ) mlix and (cosa) max are maximum values of 
H T and cos a correspondingly. 

Now, the solution of (28) is obtained in terms of Jaco- 
bian functions with modulus 



p = cos 



and period 
T 



fa 7r\ , fa 7r\ 
(s^Ms"*)' (32) 

\/2 33 K(p) 



Z{l-p + p*y ^cosd -f)V7 



(33) 
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where if (/x) is the complete elliptic integral of the first 
kind. 



We present explicit expressions for the amplitude 
squares. The convention used here is that the ampli- 
tudes are positive, which is the generic situation when 
Ht 7^ . In this convention, when Ht = the indi- 
vidual phases have discontinuities in time to account for 
the amplitudes' sign changes. The amplitude squares are 
proportional to the modes' energies and can be of great 
use for physical applications: 
2 



cm = -p 



ZT 



sn 



2K(p) 



(t-t ) 

T 



+!f (2-p + 2^1-p + p2 cos(fj) 



Cf (t) = -p 



C$(t) = p 



sn 



2K(p) 



ZT 

2 P - 1 + 2 y/1 - p + p 
2 



(t-t ) 



(34) 



ZT 



sn 



p 



(t-tg) 

T 



p 2 COS 



•/' 

(f) 



where sn(-,/i) is Jacobian elliptic function and <o is 
given in terms of the initial conditions for the ampli- 
tudes C?(0), Cf (0), C|(0) and i is defined by the initial 
conditions as: 



<o = sign(cos(^(0)) 



T 



2K(p) 



x F (arcsin y/xo, p) , (35) 



where 



x 



\/3cos(f + f) 



Z 2 T 2 (C 2 (0)-C 2 2 (0)-C 2 (0)) 
12 pK(p) 2 

(36) 

and F(-,p) is the elliptic integral of the first kind. 

Notice that each equation in (34) is a sum of two terms 
where the left terms are time-dependent and the right 
terms are not. Each right term, for instance 



T -f (2 - p + 2 ^l-p + p 2 cos v 

can be written explicitly as a function of conserved quan- 
tities -Z13, I23, Ht (expressions for p and a are given by 
(30), (31)) and is, therefore, defined by the initial condi- 
tions. 

The same is true for p and T as it follows from (32) 
and (33). In particular, one can use the equations in 
(34) to determine the minimum and maximum accessi- 
ble values of each amplitude (using the fact that sn 2 
oscillates between and 1). The characteristic energy 
variation of any resonant mode E mo d e , between these 
minimum and maximum values, has a very simple form: 
E m ode{t) ~ sn 2 (kt, p). 



D. Solution for dynamical phase 

The dynamical phase satisfies an evolution equation: 
<p = -ZH T (Ci 2 + C 2 - 2 -C 3 - 2 ). (37) 
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The solution for the dynamical phase cannot be ob- 
tained by simply replacing the solution for the amplitudes 
in the Hamiltonian Ht — C\C 2 C 3 sin tp and solving for 
tp. The reason is that non-zero tp generically evolves be- 
tween and 7r, crossing the value tp = tt/2 periodically. 
This implies that sin -1 is double- valued and thus it is 
not possible to obtain in a unique way. 

Another way to obtain the solution for dynamical 
phase might be integrating (37) in time, using the solu- 
tion for the amplitude squares, (34), but this way is also 
rather involved. On the other hand, some simple consid- 
erations allow us to find an analytical expression for the 
dynamical phase. Indeed, let us rewrite (28), taking into 
account that 4Cf = 2C\C\ and H T = C1C2C3 sintp: 



dt 



C\ = 2ZCiC 2 C 3 cos tp = 2ZH T cot tp . 



This equation can be solved for tp in each of the disjoint 
domains (0,7r) and (— 7r,0): 



tp(t) = sign(y;(0)) arccot 



sig 



n(¥>(0)) f t Cl 



2ZH T 



using the convention that the function arccot takes val- 
ues on (0, 7r). Using solution (34) together with the iden- 
tity sn'(.T, fi) = cn(x, /i)dn(x, fx) we arrive at an explicit 
expression for the dynamical phase: 



tp(t) = sign(<^o) arccot 



where 



i 1 f2K(fx)\ 



\H T \ V ZT J 



(38) 



y = sncndn 2K(fi) 



(t - t o) 
T 



,//. 



and sncn dn(-, fx) = sn(-, fx) cn(-, fx) dn(-, fx) . 

The restriction to the domain tp e (— 7r,0) U (0,7r) is 
quite general: if tp is initially in the domain (n tt, (n + 
1) 71") , n € Z, one can take tp to either (—it, 0) or (0, it) by 
an appropriate shift of 2 m tt , m G Z , without changing 
the evolution equations. Due to its special dynamics, the 
phase will remain in the domain where it was initially. 

We can simplify (38) using nome q of elliptic function 
defined as 

K'(fx) 



q = exp —tt 



Since 



K(fx) 



2K(u) / T ,, , (t-t ) 



T 



/I ^„ 1 



(2n+l)/2 



7 2n+ 



- sin(2n + l)ir 



(t - t ) 
T ■ 



(39) 



one can compute sncndn(i,/i) using the fact that 
sn' = cn dn. Explicit expression would be then a prod- 
uct of two infinite sums [25] and can be rewritten as a 



Fourier series. Since the nome q is an explicit function 
of the initial conditions, to have Fourier representations 
of our dynamical variables might be useful regarding ap- 
proximate solutions. 



E. Dynamical invariant 

Below we present a dynamical invariant for (28) which 
has been used for the constructing the solution (34). Re- 
call that a dynamical invariant depends on time, ampli- 
tudes and phases: S(t,Ci,C2,Ca,tp), with the property 
that it is a constant along any solution of the dynam- 
ical system: S(t,C 1 (t),C 2 (t),C 3 (t),tp(t)) = const. Vi. 
Generically, only local expressions can be obtained for a 
dynamical invariant, due to the multi-valuedness of the 
inverse functions involved. In this particular case, how- 
ever, since we know the period of any trajectory, this 
multi-valuedness can be eliminated partially by patching 
appropriately local expressions yielding 



S(t,C 1 ,C 2 ,C 3 ,tp) = t- 

2(t-t ) + r 

(-i)L * \t 

+ 2KiJT) 



2(t-t )+T 
2T 



T 



F (arcsin y/xt, fi) 



(40) 



where [-J is the floor function and x t can be obtained 
from the expression (36) for xq by substituting Cj(t) in- 
stead of Cj(0), for all j = 1, 2, 3. 

This dynamical invariant satisfies 

S(t,C 1 (t),C 2 (t),C 3 (t),ip(t)) = t Vi, where to is 
given in equation (35), and is an improvement of the 
corresponding formula presented in [4]. 

In Fig. 6, upper panel, we show, for fixed I\ 3 = 
2.00 and I 23 = 2.06: isosurfacc of conservation law 
Ht = 0.763, isosurfacc of dynamical invariant S = 2.69, 
solution trajectory and combined plot, in the domain 
(Ci,<p,t). 

In the middle of the upper panel, the surface S = 2.69 
is a helicoidal surface revolving around a vertical axis. 
This axis is the surface's natural interior boundary: the 
constant-in-time trajectory corresponding to the highest 
possible value of \Ht \ for given J13, 1 23 (obtained from the 
condition cosa = —1). In the present case, the highest 
possible value of \Ht\ is 1.114. This trajectory is phys- 
ically interpreted as 'maximum interference', due to the 
fact that the modes do not interact. The dynamical phase 
is constant: tp(t) = ir/2Vt, and all amplitudes are con- 
stant as well: from the condition /i = and equations 
(34), we obtain in this case: Cf (t) = 1.33, Vi . The exte- 
rior boundary of the surface is the piecewise continuous 
trajectory corresponding to the limit Ht = : in this 
limit the surface becomes non-differentiable at the 'cor- 
ners' C\ = 0, 7i3, tp = 0,7T, due to the fact that the dy- 
namical phase tp is only piecewise continuous for Ht = 0. 
This trajectory corresponds to the usual case treated in 
textbooks, when amplitudes are considered real and in- 
dividual phases vanish. 
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By looking at this figure we notice that the period 
T decreases with increasing Ht'- the trajectories closer 
to the exterior boundary are more elongated than the 
trajectories closer to the interior boundary. In fact, from 
formula (33) one can prove this property analytically. In 
Fig. 7 we plot the period T as a function of Ht- Wc 
observe in this case a reduction of the period by a factor 
0.5 when Ht is changed from to -ff max = 1.114 . 

In Fig. 6, lower panel, on the left and in the middle, 
combined plots are shown to clarify that the solution tra- 
jectory is the intersection of the isosurfaccs of Hamilto- 
nian and dynamical invariant. On the right, we show a 
combined plot of level surfaces of Manley-Rowe conser- 
vation laws /13, I 23 in the domain {C\ , Cf, Cf). 



F. Special case Ht = 

Direct substitution shows that if we put Ht = 0, then 
new modulus and period take the form 

1 rp 2K(p) 
/i = p and 1 = 

Zy h 3 

correspondingly, while the solutions for the amplitude 
squares read 

'(5?(t) = dn a ((t-to) ZVh3~,p)I 13 
< C 2 {t) = cn 2 ((t - t ) Z p) I 23 (41) 

Ci{t) = SB 2 ((t-t ) Zjh 3l p)I 23 . 

As for the dynamical phase, from Eq.(38) it is seen that in 
the limit Ht — >■ it behaves as a step function, jumping 
from to 7r sign((/?(0)): 

efl^dsyo)) (i-(-i)^jv 

To understand the meaning of this behaviour, notice 
that the Hamiltonian Ht is vanishing for ip = mr, n £ 
Z. The abrupt jumps of the dynamical phase is due to 
the jumps of the individual phases (solution not shown). 
These jumps replace the changes of sign of the modes' 
amplitudes in the usual textbook descriptions. 

As it was shown in [5], initial dynamical phase not 
in Z 7r substantially affects the magnitudes of resonantly 
interacting modes during the evolution, not only in a 
triad but also in a butterfly. This fact might have im- 
portant implications (see [5], Discussion), for instance, 
for interpreting results of numerical simulations and for 
performing laboratory experiments. 



V. GENERIC CLUSTERS 

It is well-known (see, for instance, [20, 22, 29, 30, 33? 
], etc.) that in three- wave resonance systems the most 



frequently met clusters are isolated triads or clusters con- 
sisting of two variously connected triads. Below we clas- 
sify all possible two-triad clusters - kite, butterfly and 
ray - and show how to construct new CLs making use of 
the notion of dynamical invariant. In the last subsection, 
another method is briefly outlined which was presented 
in [59, 60] and allows one to prove, in some cases, intc- 
grability of bigger clusters. 

A. Kite 

A kite consists of two triads a and b, with wave ampli- 
tudes Bj a , Bjb, j = 1,2,3, connected via two common 
modes. One can point out 4 types of kites according to 
the properties of connecting modes: PP-PP, AP-PP, AP- 
AP and AA-AA kites. In this section, PP-PP-kite with 
Bi a = B\ b and B 2a = B 2b is taken as a representative 
example. Its dynamical system reads: 

B\ a — B2 a (Z a B 3a + Z b B 3b ) , 

< B 2a = B* la (Z a B 3a + Z b B 3b ) , (42) 
B 3a = —Z a B\ a B2a , B 3 b = —ZbBi a B 2a ■ 

It has 5 conservation laws (2 linear, 2 quadratic, 1 cubic): 

{L R = Re(Z b B 3a - Z a B 3b ), Lj = lm{Z b B 3a - Z a B 3b ) , 
ha b = \B la \ 2 + \B 3a \ 2 + \B 3b \ 2 , 
hab = \B 2 a\ 2 +\B 3a \ 2 + \B 3b \\ 
Hk = lra(Z a Bi a B2 a B 3 t a + Z b B\ a Bi a B\ h )) . 

It has a dynamical invariant that is essentially the same 
as for a triad, S, after replacing 

Z = Z a + Z b , Ht = Hk{Z 2 + Z 2 )/Z 3 , 
hs = h ab (Z 2 a + Z 2 )/Z 2 - (L 2 R + Lj)/Z 2 , 
h 3 = hab{Z 2 a + Z 2 )/Z 2 - (L% + L 2 )/Z 2 . 

B. Butterfly 

A PP-butterfly consists of two triads a and b, with 
wave amplitudes Bj a , Bj b , j = 1,2,3, and connecting 
mode, say B\ a = B\ b = B\, is passive in both triads. 
The dynamical system for PP-butterfly reads 

B\ = Z a B2 a B 3a + Z b B2 b B 3b , 
B 2a = Z a B\B 3a , B 2b = Z b B*B 3b , (43) 
B 3a = —Z a B\B 2a , -83b = —Z b B\B 2b . 

We have studied this in [4] ; we present here the results 
in order to compare the dynamics of different butterfly 
types and confirm the qualitative analysis given in [30]. 
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FIG. 6: Color online. Upper panel: Triad system in coordinates (Ci , ip, t), with fixed values of the Manley-Rowe conservation 
laws: Ji3 = 2.00 and I23 = 2.06. Left: level surface of conservation law Ht = 0.763, Eq.(29). Middle: level surface 
of dynamical invariant S = 2.69, Eq.(40). Right: solution trajectory (Ci(t),ip(t)), Eqs.(34),(38), corresponding to Ht = 
0.763, S — 2.69. Lower panel: Left: combined plot of level surface of dynamical invariant S = 2.69 and solution trajectory 
(Ci(i) 2 , <p(t)). Middle: combined plot of level surface of conservation law Ht = 0.763, level surface of dynamical invariant 
S = 2.69 and solution trajectory (Cj (t), </?(£))■ Notice the general property that the intersection of the level surfaces of Ht and 
S is the solution trajectory. Right: combined plot of Manley-Rowe conservation laws /13 = 2.00 and I23 = 2.06 , in coordinates 

(cicid). 



Sys.(43) has 3 quadratic CLs analogous to (5) and 1 cubic 
CL corresponding to its Hamiltonian: 



>3a| 



^236 = I-B2&I + 
! 3b| 2 I 

Hpp = lm(Z a B 1 B2 a B^ a + ZbB 1 B 2 bB^ b ) 



3361 



I ab = |Bl| 2 + |i?3a| 2 



(44) 



Standard amplitude-phase representation. Here, one 
can rewrite the cubic conservation law as 



Hpp — Cpp (Z a C2aC^ a sin ip a + ZbC 2 bC 3 b sin ipb) . 



Here 



(45) 



fa = Ola + ' 



^3a, fb 



lib 



12b 



are dynamical phases and Crpp) is the real amplitude 
of a common mode in PP-buttcrfiy. This allows us to 



reduce Sys.(43) to only four real equations: 



C*3q = — ZaPppCia COS <fi a , 

C 3 b = —Zb CppC 2 b cos tpb , 



ip a = -Hp P (C 

fib = -H PP (C 



pp 
pp 



Co 



c 3 -a% 



(46) 



s < — 1 s~y — 2 

°2b _ °36 



Now the overall dynamics of the PP-butterfly is confined 
to a 3-dimensional manifold. The same can be done for 
the two other types of butterflies. 

Modified amplitude-phase representation. The follow- 
ing change of variables was suggested in [4] : 

a a = arctan(C 3a /C2a) , «b = arctan (C 3 b/C 2 b) ■ 
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T(H T ) 



T\ Ht =o =3.51 x 10' 




3.7 



7.4 
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FIG. 7: Color online. With fixed values of the Manley-Rowe 
conservation laws: /13 = 2.00 and 723 = 2.06 , plot of period T 
from Eq.(33) as a function of Ht, normalized with respect to 
the period at Ht = 0. The decreasing character is a generic 
feature of this function. 
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FIG. 9: NR-diagram for PP-butterfly. 



independent CLs will be suitable. For instance, in the 
case Hpp = we could choose the conserved quantity 

A5 = sin tpi, sin 2a& instead of A a = sin ip a sin 2a a 

(see Table I, 1.2 PP) but not both because they are func- 
tionally dependent: 



Z a l23aA a + ZhlllhA 



0. 



Generally, we try to find the simplest presentation for 
our new constants of motion. 



FIG. 8: NR-diagram for PP-kite. 



with the inverse transformation being 

I C2a = Vh3a COS CL a , C^ a 

I C 2 b = \fhzb cos a b , C 3b = 



(47) 



This change of variables allows further substantial sim- 
plification of (43) and (46): 



C. Ray 

Analogously with the previous case, AA-butterfly is 
a two-triad cluster with a common mode which is A- 
mode in both triads, B 3a = B 3 b- Dynamical system and 
Manley-Rowe constants read: 




Z a B2 a B 3a , Bib — ~Zt,B2 b B 3a 

^a^la^3a , &2b — ^h£> lf) ±>3 a 
— Z a Bi a B2a — ZbBn,B2b ■ 



I\2b 



\B U 



\B 



2h\ 



(51) 



(52) 



a a = -Z a C PP cos ip a , 
a b = -Z b Cpp cos ipb , 

<p a = Z a Cp P (cota a - tana Q )sin^ a - H PP /Cpp, 
(fib = ZbCpp (cot ab — tancvh) sin^b — Hpp/C PP . 

In these new variables, the amplitude Cpp > reads 



(48) 



Cpp = y I a b- h3a sin 2 a a - h3b sin 2 a b 
and the Hamiltonian is now 
Cpp 



(49) 



Hpp — 



-Z a I 2 3a sin ip a sin2a a 



C 



pp 



Z 6 / 2 36 sin (fib s'm2a b . 



(50) 



Eqs. (48)-(50) represent the final form of our three- 
dimensional general system in the modified amplitude- 
phase presentation. 

A few cases of integrability in quadratures of the PP- 
butterfly were presented in [4]. The results are collected 
in Table I below. Of course, the form of the conservation 
laws is arbitrary in the sense that any set of functionally 



FIG. 10: NR-diagram for AA-butterfly. 

The integrability of (51) can be investigated along the 
same lines as for (43) above. The analysis is omitted 
here. We just partly outline one particular case of this 
cluster: AA-r&y, which can be regarded as a degenerate 
AA-butterfly, so that = u)2b = ^3/2. 

In this case, the dynamical system obtained from first 
principles will have the form 



B\ a — ZaB^aB-i , Bb 



ZbB^B 3 



B2a — Z a B* la B 3 , i? 3 — —Z a B\ a B2a — 2ZbB b 



(53) 



Notice that there is a factor 2 in the last term of last 
equation, which would not appear if we made the direct 
substitution B^ = B 2 b = Bb into system (51). Rather, 
the simple change of variables Bib = B 2 b = V^B b will 
transform the AA-butterfly (51) into the ray equations 
(53). This means in particular that integrable cases of 
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AA-buttcrfly can be directly mapped to some integrable 
cases of AA-ray. Another interesting point is that AA- 
ray cluster might also have a nice mechanical model - 
Wilberforce pendulum [38], the problem is presently un- 
der the study. 

Conservation laws for AA-ray are inherited from con- 
servation laws for AA-buttcrfly: 



I ab = \B la \ 2 + 2\B b \2 + \B 3 \ 2 



(54) 



I\2a = \B\a\ 
Lib = \Bla\ 2 

Hra,y = lm(—Z a B\ a B2 a Bz — 2Z b B\B 3 ), 
with dynamical phases 

tfia = Q\a + #2a — #3, = 29 b — #3. 



This reduces four complex equations (53) to only four 
real ones: 



Cla = Z a C 2a C 3 COS (f a , 

C b = Z b C b C 3 cos tp b , 

(p a = -Z a C 3 (g^ + g^) sinp„ + H lay /C\ 



(55) 



^ Lp b = -2Z b C 3 sin ip b + H ra . y /C§ , 
with Hamiltonian 

H lay = -C 3 (Z a Ci a C 2 a sin ip a + 2Z h C\ sin ip b ) 



(56) 



in terms of the amplitudes and phases. 

Consider the simple case when initially ip a = fb = 0. 
Then H lay = 0, phases remain zero for all times, and the 
equations of motion reduce to 



Cla = Z a C2aC 3 , 
C*2a = Z a C\ a C 3 , 

Cb = Z b C b C 3 , 

C 3 = —Z a G\ a C2a — ZZbCfr 



with two Manley-Rowc constants of motion 
I I\2a = C la - C| a , 

I I ah = C{ a + 2C\ + Cf . 



(57) 



and a new one 



2Z a \nC b + Z b ln 



C 2a — C\a 
C 2a + Cla 



(58) 



(59) 



where notations 

Cia = \J l\2a cosh a , C 2a - 
are used. 

Notice that this way only non-polynomial additional 
conservation laws are obtained (see Table I). 




/ \ A 

J, 



FIG. 11: NR-diagrams for 3-star clusters. From right to left, 
from up to down: 3-star- A (three A-connections) , 3-star-P 
(three P-connections) , 3-star- 1A-2P (one A-connection and 
two P-connections), 3-star-2A-lP (two A-connections and one 
P-connection) 



D. Star 

A cluster of N triads, all connected via one common 
mode is called N-stai cluster. Again, intcgrability of TV- 
star depends on the types these connecting modes have 
in each triad of a cluster. NR-diagrams for all possible 
types of 3-stars are shown in Fig. 11. 

iV-star cluster is probably the only presently known 
type of cluster for which an analytical study has been 
performed for arbitrary finite number N. The main idea 
can be briefly formulated as follows. ./V-star cluster has 
2N+1 degrees of freedom, N+l Manlcy Rowc constants 
of motion and one Hamiltonian, that is, we already have 
N+2 independent first integrals in involution. To find 

— 1 additional integrals of motion, one can use con- 
struction of Lax operators, Painleve analysis and irre- 
ducible forms (see [45, 46, 59, 60], etc.; terminology used 
therein is pump and daughter wave for A- and P-mode 
correspondingly) . The dynamical system, say for JV-star- 
A, is regarded in the form 

Bij = i X 3 B 3 B 2j , i! 2 , = i A, /.',/>; ,. (60) 

N 

B 3 = iJ2^BijB 2j . (61) 

j=i 

Additional conservation laws found this way have neces- 
sarily polynomial form. The results for a generic N-star 
cluster are as follows: N-star-A (with all A-connections) 
and N-star-P (with all P-connections) are integrable for 
arbitrary initial conditions if 



1 or 2, 



examples of corresponding NR-diagrams shown in the 
Fig. 11, upper panel, for N = 3. A-star cluster with 
mixed A- and P-connections (Fig.ll, lower panel) has no 
additional polynomial conservation laws. Complete set 
of additional polynomial conservation laws for integrable 
N-star cluster is omitted here for sake of place, and it can 
be found in [60]. Example for the case of AA-buttcrfly 
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Conditions 




Additional CLs 


fa = fb = 


1.1. 

PP-but. 


Z b arctan (^j^ j — Z a arctan (j^J 


(fa, <fb 7^ 0, ffpp = 


1.2. 
PP-but. 


A a = sin y> a sin 2a a , with a a = arctan (C3a/C2a) 




1.6. 
PP-but. 


I 1 + z °)arccos 1 1 (1 + z " a ) arccos 1 ^—^ 1 
with at = arctan (C3b/C2b) and A;, = sin </P{, sin 2ab 


ffpp 7^ 0, Z a = Zj, 


1.4. 
PP-but. 


ClaCla + C\ b C\ b + 2C2aC3aC 2b C 3b COs((fi a - lf b ) - C?(C? - C"| a + C| a - C\ b + C| 6 ) 


= <Pb = 


2.1. 

AA-ray 


2Z a In a + Z b In [(C 2a - C la )/(C 2a + C la )] 

With Cla = V^12a COsh Q , C 2 a = V^12a SUlh Q , /l 2 a = C\ a ~ C| a 



TABLE I: Examples of non-polynomial conservation laws. 



with Z a — 2 reads 

AiBxBrflBt + BlBlBiB^BxBl + B 2 B* 2 ) 
-2{B 3 B{B* 2 + BtBtB 2 ) 2 - + B 2 B*) 2 

+AB 1 B* 1 B 2 B* 2 ]{B i Bl + B 5 B*). (62) 

The general way to investigate intcgrability of a generic 
cluster would be to apply the theory of normal forms (see, 
for instance, [50, 52]) for each dynamical system describ- 
ing a resonance cluster is a normal form. But of course, 
most generic clusters demonstrate chaotic behavior and 
numerical investigations are unavoidable. 

Here we outline briefly which facts are important in 
order to perform sensible numerical simulations with res- 
onance clusters. The fact that our systems are Hamilto- 
nian, allows us to perform numerical simulations based 
on the Hamiltonian expansion of the corresponding dy- 
namical system. Poincare sections are well-known in- 
struments to show clearly whether or not a dynamical 
system demonstrates chaotic behavior, so all the results 
of our numerical simulations arc illustrated by the cor- 
responding Poincare maps. One of the most tedious and 
time-consuming parts of the corresponding simulations is 
the choice of initial conditions. A special procedure has 
been worked out [43] that guarantees a uniform distribu- 
tion of initial conditions according to Liouville measure, 
and assures as well that all conservation laws have the 
same value on each Poincare section. In Fig. 12 an ex- 
ample of Poincare section is shown. To trace effects due 
to the dynamical phases one has to use amplitude-phase 
representations, of course. 




0.5 1 1.5 2 2.5 3 



FIG. 12: Example of Poincare section for PP-butterfiy with 
Z a /Z b = 3/4. 



VI. COUPLING COEFFICIENT 

As we have shown above, the integr ability of a res- 
onance cluster depends on the magnitude of the corre- 
sponding coupling coefficients. Expressions for coupling 
coefficients in canonical variables has been deduced for 
various types of wave systems possessing three- wave reso- 
nances: irrotational capillary waves [66]; rotational capil- 
lary waves [10]; irrotational gravity-capillary waves [47]; 
rotational gravity-capillary waves [18], drift waves [56], 
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etc. They usually have a nice compact form, for instance, 
coupling coefficient for irrotational gravity-capillary 
water waves reads 

(w|-w 2 w 3 +w 2 ) 

fci - uj 2 k 2 + w 3 fc 3 , 

where w = (.gfc + erfc 3 ) 1 / 2 and and c are gravity ac- 
celeration and surface tension correspondingly. However, 
transformation of these expressions from the canonical to 
physical variables is not an easy task. 

On the other hand, the application of multi-scale meth- 
ods yields expressions for the coupling coefficients di- 
rectly in physical variables. For instance, coupling co- 
efficients of the system of three resonantly interacting 
atmospheric planetary waves, with ui ~ m/[n{n + 1)], 
have the form [33] 

Z[n 2 (n 2 + 1) - n 3 (n 3 + l)]/[ni(m + 1)], 
Z[n 3 (n 3 + 1) - m(m + 1)]/M"2 + 1)], (63) 
Z[n 2 (n 2 + 1) - m(m + l)]/[n 3 (n 3 + 1)], 

with 

y-Tr/2 "V dtp dip 

Here two spherical space variables are the latitude tp, 
— ir/2 < tp < 7r/2, and the longitude A, < A < 2tt, 
and the notation P^ is used for P^ 3 (sin tp) which is the 
associated Legendre function of degree rij and order rrij. 



The multi-scale method is quite straightforward and 
can be programmed in some symbolical language [33]. 
However, only numerical magnitudes of the coupling co- 
efficients have been computed for selected solutions of the 
resonance conditions, not an explicit algebraic formulas. 
The problem is due to some "bags" in Mathematica in 
computing integrals of the form 

I sin(mx) sm.(nx)dx, with to,tt,gN, (64) 
Jo 

more discussion can be found in [33]. 

For completeness of presentation we show below a part 
of the formula for one coupling coefficient, obtained by 
a combination of symbolical programming in Mathemat- 
ica and automatic search-and-replace of formulas of the 
type (64) - which Mathematica cannot handle directly - 
for ocean planetary waves with vanishing boundary con- 
ditions, w ~ (m 2 + n 2 ) -1 / 2 , in physical variables. The 
general form of the coupling coefficients is 

a,/[(TO 2 + n 2 K 2 ], J = 1,2,3, (65) 

with coefficients a.j being functions of the wavenumbers 
mi, rii, ■ ■ ■ , m 3 , n 3 . Notations E and I below are for ex- 
ponent exp and imaginary unit i, all other notations are 
self-explanatory. Notation (...) at the end is used for 
eight more pages, needed to accomplish this formula. 



(E~(-I (ml + m2 + m3 + Sqrt[ml"2 + nl"2] + Sqrt[m2~2 + n2~2] - Sqrt[m3"2 +n3"2] )\ [Pi] ) 
((I (m2~2 m3 n2 - 2 m3~3 n2 + m3 n2"3 + m2 m3 n2 Sqrt[m2~2 + n2~2] + 2 m2~3 n3 - m2 m3~2 n3 
+ 2 m2 n2~2 n3 + 2 m2~2 Sqrt[m2~2 + n2~2] n3 - m3~2 Sqrt[m2~2 + n2~2] n3 + n2~2 Sqrt[m2~2 
+ n2~2] n3 - 2 m3 n2 n3~2 - m2 n3~3 - Sqrt[m2~2 + n2~2] n3~3 + m2~2 n2 Sqrt[m3~2 + n3~2] 

- 2 m3~2 n2 Sqrt[m3~2 + n3~2] + n2~3 Sqrt[m3~2 + n3~2] + m2 n2 Sqrt[m2~2 + n2~2] Sqrt[m3~2 
+ n3~2] - m2 m3 n3 Sqrt[m3~2 + n3~2] - m3 Sqrt[m2~2 + n2~2] n3 Sqrt[m3~2 + n3~2] 

- n2 n3~2 Sqrt[m3~2 + n3~2]))/(ml + m2 - m3 + Sqrt[ml~2 + nl"2] + Sqrt[m2~2 + n2~2] 

- Sqrt[m3~2 + n3~2] ) + (I (m2~2 i3 n2 - 2 m3~3 n2 + m3 n2~3 + m2 m3 n2 Sqrt[m2~2 + n2~2] 
+ 2 m2~3 n3 - m2 m3~2 n3 + 2 m2 n2~2 n3 + 2 m2~2 Sqrt[m2~2 + n2~2] n3 - m3~2 Sqrt[m2~2 

+ n2~2] n3 + n2~2 Sqrt[m2~2 + n2~2] n3 - 2 m3 n2 n3~2 - m2 n3~3 - Sqrt[m2~2 + n2~2] n3~3 
+ m2~2 n2 Sqrt[m3~2 + n3~2] - 2 m3~2 n2 Sqrt[m3~2 + n3~2] + n2~3 Sqrt[m3~2 + n3~2] 
+ m2 n2 Sqrt[m2~2 + n2~2] Sqrt[m3~2 + n3~2] - m2 m3 n3 Sqrt[m3~2 + n3~2] - m3 Sqrt[m2~2 
+ n2~2] n3 Sqrt[m3~2 + n3~2] - n2 n3~2 Sqrt[m3~2 + n3~2]))/(ml - m2 + m3 - Sqrt[ml~2 
+ nl~2] - Sqrt[m2"2 + n2~2] + Sqrt[m3~2 + n3~2]) - (I (m2~2 i3 n2 - 2 m3~3 n2 + m3 n2~3 

- (...) 



VII. SUMMARY 

(..to be added: numerical simulations, NR-reduced 
models..) 
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